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20.  ABSTRACT  (continued) 


sinusoidal  amplitude).  / While  the  deterministic  approach  is  valid  only  when  the 
initial  conditions  are  known  exactly,  it  must  be  superseded  by  a stochastic 
model  when  these  conditions  are  known  only  in  a statistical  sense.  This  model 
is  successfully  applied  to  the  breakup  of  liquid  jets  where  turbulence  is 
assumed  to  provide  the  random  disturbance. 
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I.  INTRODUCTION 


In  classical  stability  analysis,  one  observes  the  behavior  of 
small  departures  from  the  equilibrium  state.  These  disturbances  decay 
in  stable  systems,  retain  their  initial  amplitude  in  neutrally  stable 
systems,  and  are  amplified  in  the  unstable  case.  It  is  often  conven- 
ient to  choose  these  small  perturbations  from  a set  of  orthogonal  and 
complete  functions  since  in  the  linear  case  any  arbitrary  initial  per- 
turbation can  be  expressed  as  a linear  combination  of  members  from  this 
set  so  that  one  may,  without  loss  of  generality,  confine  one's  analysis 
to  an  initial  perturbation  which  is  described  by  a single  but  arbitrary 
member  of  this  set.  In  many  cases,  the  set  par  excellence  consists  of 
sinusoidal  functions.  In  addition  to  their  obvious  mathematical  con- 
venience, they  are  easily  reproduced  experimentally  so  that  a direct 
comparison  between  theory  and  experiment  is  possible.  In  these  experi- 
ments, the  initial  perturbation  is  usually  too  small  to  be  directly  per- 
ceptible. As  it  grows  in  time,  it  manifests  itself  macroscopically,  the 
effect  of  the  instability  being  similar  to  an  amplifier  for  this  small 
disturbance.  Such  controlled  experiments  are  characterized  by  a high 
degree  of  regularity  which  disappears  as  soon  as  the  disturbance- 
producing  mechanism  is  removed.  Here,  the  unstable  system  obtains  its 
initial  perturbation  from  the  random  ambient  vibrations  always  present 
in  its  environment.  Consequently,  the  macroscopic  behavior  of  this  sys- 
tem is  also  random.  A deterministic  description  of  this  process  is 
ruled  out  since  the  initial  perturbation  is  not  known  exactly.  However, 
a knowledge  of  the  statistical  and  spectral  properties  of  the  initial 
perturbation  is  sufficient  to  determine  the  statistical  and  spectral 
properties  of  the  unstable  system  after  the  initial  perturbation  has 
been  applied. 

This  paper  discusses  the  response  of  an  unstable  linear  system  to 
stationary  gaussian  noise  with  zero  mean.  The  treatment  is  very  similar 
to  that  of  calculating  the  output  of  a linear  filter  with  stochastic 
input  but  with  one  major  difference:  the  present  analysis  involves  both 
space  and  time  variables.  The  random  functions  describing  system  param- 
eters are  stationary  and  gaussian  in  the  space  coordinates,  but  generally 
increasing  functions  of  time  (because  of  the  instability).  We  show  that 
when  a certain  parameter  (say,  a flow  velocity)  is  perturbed  with  sta- 
tionary gaussian  noise  with  a given  spectral  density,  the  perturbation 
grows  in  time  but  remains  gaussian  - its  standard  deviation  being  a func- 
tion of  time.  Furthermore,  when  the  spectral  density  becomes  confined  to 
a narrow  band  of  wave  numbers,  the  amplitude  of  the  perturbation  is 
Rayleigh  distributed.  The  present  model  is  found  to  be  in  good  agreement 
with  the  data  of  Phinney1,  suggesting  a simple  model  for  the  breakup  of 
turbulent  liquid  jets;  one  may  view  the  turbulence  as  a source  of  random 


* 


1. 


R.  E.  Phinney,  " The  Breakup  of  a Turbulent  Liquid  Jet  in  a Gaseous 
Atmosphere,"  J.  Fluid  Meah. , Vol.  60,  1973,  p.  689. 


perturbations  for  the  liquid  jet.  Since  the  turbulence  amplitude  is 
small,  one  could  expect  the  linear  breakup  theory  of  Rayleigh2  to  be 
applicable,  with  modifications  described  in  the  present  report  to  take 
into  account  the  stochastic  nature  of  the  perturbation. 


II.  LINEAR  THEORY 


In  linear  stability  theory,  one  can  often  decompose  perturbations 
of  the  field  variables  (e.g.,  velocity  components,  pressure,  movable 
boundary  coordinates)  into  elementary  components  of  the  form 

exp  i(kx  - uit) 

where  x and  t are  space  and  time  coordinates  respectively,  k = 2ir/X  is 
the  wave  number  of  that  component  (with  wavelength  X)  of  the  perturba- 
tion applied  at  time  t = 0,  and  w(k)  = Wj(k)  + iw.,(k)  is  the  complex 

frequency.  We  can  express  a typical  field  variable,  denoted  by  tl<  (x,t) 
by  a linear  superposition  of  these  elementary  components: 


ij/(x,t)  = /k  H'(k) 


Hkx-u.jt)  u t dk 
e e 


(1) 


where  vCk)  is  the  Fourier  transform  of  the  initial  perturbation.  Eq.  (1) 
can  be  written  in  the  form 


«Kx,t)  = /k  F(k,t)  S'(k)  elkx  dk 


(2) 


where 


-iu>  (k)t  u>,(k)t 

F(k,t)  = e 1 e (3) 

This  form  is  reminiscent  of  a linear  filter  with  input  waveform 
ip (x , o)  at  t = 0 and  output  waveform  iHx.t)  at  time  t.  F(k,t)  is  the 
time-evolving  Fourier  transform  of  the  impulse  response. 

When  the  input  (initial  condition)  is  stochastic,  then  the  output 
<J>(x,t)  will  also  fluctuate  randomly.  For  the  sake  of  definiteness,  and 
mathematical  tractability,  we  consider  the  initial  conditions  to  have 
the  following  properties  in  the  stated  coordinate: 

a.  Stationary 

b.  Gaussian  with  zero  mean 


2. 


Lord  Rayleigh , Theory  of  Sound,  Vol.  II,  Dover,  New  ,'ork  (1945) 
pp • 351—355. 
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c.  Ergodic 

The  system  is  assumed  to  be  linear.  Time  plays  the  role  of  a family 
parameter  which  we  consider  fixed  for  any  input  - output  linear  oper- 
ation. Since  the  system  is  linear,  a gaussian  input  will  result  in  a 
gaussian  output.  We  characterize  our  random  process  by  its  autocorre- 
lation function.  This  characterization  represents  a simple  statistical 
average  based  on  the  second  order  probability  distribution  function  for 
the  process.  The  autocorrelation  function  is  defined  by 


R(y*t)  = lim 


iHx,t)  <Hx+y,t)dx 


(4) 


# z-**»  -z 

The  spectral  density  for  this  process  is  obtained  from  the  Wiener- 
Khintchine  relation 


G(k,t  ) = 2^/  R(y,t) 


»-iky 


dy 


whence 


R(y,t)  = / G(k,t)  elky  dk 


(5) 

(6) 


We  now  write  the  input  - output  relation  for  the  spectral  densities 

G(k,t)  = |F(k,t)|2  G(k,o)  (7) 

where  the  input  is  applied  as  an  initial  condition  with  spectral  density 
G(k,0)  at  time  t = 0 and  the  output  is  expressed  in  terms  of  either  its 
autocorrelation  function  R(y,t)  or  of  its  spectral  density  G(k,t). 

Since  G(k,0)  is  given,  we  can  view  |F(k,t)p  as  a time  evolution  operator 
which  relates  the  state  of  the  system  at  time  t to  the  given  initial 
conditions.  We  can  easily  find  an  explicit  form  for  the  time  evolution 
operator.  Using  Eqs.  (3)  and  (7),  we  have 

2u>-(k)t 

G(k,t)  = e G (k, 0)  (8) 

so  that  the  spectral  evolution  depends  only  on  the  imaginary  part  of  the 
complex  frequency. 

Consider  now  the  set  of  stationary  random  gaussian  field  variables 
with  zero  mean 

<♦„>  - <9> 


3. 


W.  B.  Davenport,  Jr.,  and  W. 
of  Random  Signals  and  Noise, 


L.  Root,  An  Introduction  to  the  Theory 
McGrau  Hill,  New  York  (1958). 
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The  probability  of  measuring  ^ at  x = x^,  ^ at  x = x^>  . . . , ^ at 
x = x^,  all  at  time  t is3 


...  dtf/^  _ j 

p{*N)  d*j  ...  d*N  • ^)N/2  i,  exp  2JIJ 


N 

E I A I if/  if/ 

1 'm,n  m n 

n,m 
= 1 


(10) 


where  A is  the  co-variance  matrix  with  elements 


X = R(x  - x ,t) 
m.n  m n 


(11) 


and  co-factors  I A I . For  a stationary  process,  let  y = x - x ; also 
1 'm,n  m n 

define  o2(t)  = o2  = R(0,t).  We  then  have 


°2 (t)  R(y,t)  — R(y,t) 

R(y,t)  a2(t)  ....  R(y,t) 


A = 


R(y,t)  R(y,t) 


with 


a2(t) 


(12) 


00  2w,)(k)t  . , 

R(y,t)  = /_„  G(k,o)  e e1  y dk 


(13) 


and 


00  2u>_(k)t 

02(t)  = / G(k,o)  e dk 


(14) 


We  have  taken  advantage  of  the  fact  that  for  a stationary  process  the 
multivariate  gaussian  distribution  with  zero  mean  is  totally  specified 
by  the  auto-correlation  function,  which  in  turn  is  related  to  the  ini- 
tial condition  by  Eq.  (13). 
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III.  WIDEBAND  INPUT  WITH  DOMINANT  MODE 


We  proceed  by  studying  a particularly  simple  but  yet  often  encountered 
form  for  G(k,0)  and  u>-(k).  Imagine  the  initial  condition  to  be  gaussian 
white  noise 


G(k,0)  * 1 

and  ^2 (k)  to  be  strongly  peaked  at  k = kQ 
u>2(k)  > 0 
dw- 

-rr—  =0  at  k = k 
dk  o 

d2u>_ 

-dF<0atksko 


(15) 

(16a) 

(16b) 

(16c) 


Note  that  Eq.  (15)  represents  an  assumption  about  the  environment  while 
Eq.  (16)  is  based  upon  a description  of  the  physical  system.  In  prac- 
tice, these  quantities  can  be  measured  and  the  exact  techniques  of  the 
previous  section  may  be  applied.  In  a large  number  of  cases,  however, 
Eqs.  (15  - 16)  represent  a sufficiently  good  approximation  to  reality, 
for  example,  a liquid  jet  subjected  to  wideband  vibrations  (which  are 
assumed  to  have  Gaussian  statistics).  We  start  by  expanding  w7(k) 
about  kQ  and  keeping  only  the  first  two  non-vanishing  terms: 


u>2  (k)  = “2^0^  + 55  ' kQl  “2"(ko) 


(17) 


since  (^'(k^  = ^ can  usec*  to  calculate  the  spectral  density, 

which  takes  the  form  of  a Gaussian  distribution 

k-k 


G(k,t)  = A (t)  exp  [-h  (- 


IN  «“> 

■V] 


(18) 


with  amplitude 

2u)_(k  )t 

*Ct)  . . 2 ° 

and  standard  deviation 

Ss  = [ 2 |^2"(ko)|t]-Js 

The  auto-correlation  function  is  readily  obtained: 


(19) 


(20) 





I m ' '■ — '■  ■ 


.... 


R(y,t)  = /2n  A(t)  cos  (kQy)  exp  [ -*5(ySs)  ] 


2u.(k  ) t r n .4 

o2(t)  = e ° |w,"(k  ) 


” “ I "2  '-  QJ  1 

The  singularities  in  Eqs.  (20)  and  (22)  need  not  be  alarming;  they  merely 
indicate  that  the  spectral  and  statistical  distributions  of  the  perturba- 
tion tend  toward  uniform  distributions  as  t becomes  small.  Examination 
of  Eq.  (22)  shows  that  o (t)  decreases  from  a very  large  value  at  t = 0 
to  a minimum  at  t = l/4a)2(kQ)  = t and  thereafter  increases  without 

bound.  Typically,  x is  a very  small  quantity.  For  example,  in  a liquid 
jet,  it  is  of  the  order  of  one-fiftieth  of  the  breakup  time.  This 
singular  behavior  at  small  times  is  but  a reminder  that  white  noise  with 
infinite  bandwidth  is  unphysical. 

Note  from  Eq.  (20)  that  as  time  increases  the  harmonic  contents  of 

the  perturbation  tend  to  become  concentrated  in  a narrow  band  of  wave- 

numbers  centered  about  k = k . Functions  such  as  these,  which  have  a 

o 

sharply  peaked  spectral  density,  can  be  written  in  the  narrowband 
form1*  » 5 

>Kx,t)  = £(x,t)  cos  [kQx  - oijt  + <Kx,t)]  (23) 

where  S(x,t)  and  4> (x , t ) represent  the  envelope  and  phase  of  ij>(x,t) 
respectively,  and  are  slowly  varying  functions  of  x.  ij;(x,t)  owes  its 
almost  sinusoidal  appearance  to  the  predominance  of  spectral  components 
at  k = kQ.  The  presence  of  adjoining  components  gives  rise  to  slow 

changes  in  the  envelope  and  phase.  If  the  phase  is  uniformly  distrib- 
uted (random  phase  approximation),  then  it  can  be  shown5  that  the  fol- 
lowing statistical  distributions  hold 


= 2? 


(24a) 


f5(0  = ^7-exp  [-^(f) 2] 


(24b) 


4.  J.  S.  Bendat,  Principles  and  Applications  of  Random  Noise  Theory, 
John  Wiley  <S  Sons , New  York  ( 1958) . 

5.  M.  Schwartz,  Information  Transmission , Modulation  and  Noise , 
McGrow  Hill,  Hew  York  (1970). 


10 


V*)  - ~~  exp  [->s(£)2]  . (24c) 

v a/ 2rr 

Although  we  have  shown  narrowband  signals  to  arise  from  the  preferential 
amplification  of  spectral  components  near  the  wave  number  kQ,  it  is  easy 

to  see  that  they  can  also  be  produced  when  G(k,0)  is  sharply  peaked  at 
k = kQ,  say. 


IV.  APPLICATION  TO  TURBULENT  JET  BREAKUP 


Phinney1  has  measured  the  statistical  distribution  of  the  breakup 
length  of  turbulent  liquid  jets.  The  breakup  length  is  defined  as  the 
length  of  the  continuous  portion  of  the  jet,  measured  from  the  nozzle  to 
the  breakup  point  where  drop  formation  occurs.  For  artificially  per- 
turbed laminar  jets  (e.g.,  with  an  imposed  sinusoidal  vibration  of 
constant  amplitude) , the  breakup  length  is  a constant  whose  magnitude 
depends  on  the  amplitude  and  wavelength  of  the  imposed  perturbation  as 
well  as  on  hydrodynamic  parameters  such  as  jet  diameter,  surface  tension 
and  fluid  density.  When  no  artificial  disturbance  is  used,  or  when  the 
jet  is  turbulent,  the  breakup  length  fluctuates  about  an  average  value. 
Phinney  states  that  the  distribution  of  the  turbulent  breakup  length  is 
described  well  by  a Gaussian  but  remarks  that  there  exists  no  theoret- 
ical background  in  which  to  fit  the  observation.  We  do  not  propose  such 
a theory  here,  but  have  posed  the  following  question.  Can  we  achieve  a 
reasonably  good  correlation  with  experimental  data  by  assuming  that 
turbulent  jet  breakup  occurs  as  a result  of  random  perturbations  acting 
on  an  otherwise  laminar  capillary  jet?  We  shall  assume  that  Rayleigh's 
linear  theory2  for  the  breakup  of  inviscid  capillary  jet  applies,  with 
the  sole  modification  that  the  small  initial  perturbation  is  allowed  to 
fluctuate  randomly.  Thus,  the  role  of  turbulence  is  limited  to  that  of 
a driving  mechanism  which  initially  deforms  the  jet  surface  in  a random 
manner.  This  ruffled  surface  initiates  the  process  of  capillary  insta- 
bility which  ultimately  results  in  jet  breakup.  Rayleigh's  criterion 
for  breakup  can  be  written  as 


h = Co  exp  [w  (k) 


(25) 


where  Co  is  the  initial  perturbation  amplitude,  N = pd0U02/T  is  the 
Weber  number,  d0  is  the  equilibrium  jet  diameter,  Uq  is  the  jet  veloc- 
ity, T is  the  surface  tension  and  p is  the  density  of  the  fluid.  Eq. 
(25)  states  that  the  breakup  occurs  when  the  initial  amplitude  Co 

(rendered  dimensionless  with  vespect  to  do)  has  been  exponentially 


11 


■■■■  totmu 


amplified  to  the  size  of  the  jet  diameter.  The  breakup  time  and  jet 
velocity  have  been  combined  to  give  L,  the  breakup  length  measured  in 
jet  diameters.  The  function  w(k)  is  the  amplification  growth  rate  and 
depends  on  the  wavenumber  of  the  perturbation 


w(k)  = [ 


8T 

Pdo3 


k 


(1-k2) 


VIO  , 
I0OOJ 


(26) 


where  k is  the  dimensionless  wavenumber  and  Iq  and  1^  are  modified 

Bessel  functions.  Eq.  (25)  provides  a deterministic  relationship 
between  the  initial  perturbation  amplitude  and  jet  breakup  length.  If 
however  the  initial  perturbation  amplitude  is  not  fixed,  but  is  dis- 
tributed according  to  some  statistical  law,  say  f (C0),  then  the  break- 

Co 

up  length  will  be  distributed  according  to  fL(L),  the  two  distribution 

variables  being  related  by  Eq.  (25) . One  can  express  one  statistical 
distribution  in  terms  of  the  other  by  a simple  variable  transformation 
which  conserves  probability 


fL(L)  = f 


Co 


(27) 


Our  problem  then  is  reduced  to  that  of  finding  f (Co)-  To  this  end 

we  note  that  w(k)  has  a peak  at  k0  = 0.698.  The  effect  of  this  is  to 

preferentially  filter  initial  excitations  so  that  only  a narrow  range  of 
wavenumbers  centered  about  k0  can  significantly  participate  in  the 

breakup  process.  It  has  been  shown4’5  that  the  statistical  distribution 
of  the  amplitude  in  a narrowband  process  is 

fCo  <So)  ■ § exp  [ <5?2  1 (W 

where  o0  is  the  r.m.s.  value  of  the  random  amplitude.  Eqs.  (."’5),  (27) 
and  (28)  can  be  used  to  obtain 


fLa) 


n 

4^? 


-2QL 

e 


exp 


[- 


where 


0 = 


(£-)%  u(k0) 

w 


Eq.  (29) 


has  a maximum  at 


(29) 


(30) 
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which  is  the  mode  of  that  distribution  and  is  used  in  defining  a con- 
veniently normalized  breakup  length 


Phinney's  measured  breakup  lengths  were  normalized  to  the  mean  breakup 
length  and  were  plotted1  to  show  that  these  data  exhibit  gaussian  scat- 
tering about  the  mean.  For  comparison  with  the  present  theory,  Phinney's 
data  have  been  replotted  along  with  our  theoretical  breakup  length 
distribution.  The  result  is  shown  in  Fig.  (1),  where  the  only  free 
parameter,  o0,  has  been  adjusted  to  minimize  the  mean  square  error,  oq 

was  thus  determined  to  have  the  value  0.008,  a reasonable  figure  for  the 
r.m.s.  amplitude  of  turbulent  fluctuations;  unfortunately,  no  experimen- 
tal value  was  reported.  The  fit  is  not  uniformly  good,  but  a chi-square 
analysis  revealed  that  the  level  of  confidence  is  better  than  97°<.  for 
0.85  < t < 1.15  and  drops  rapidly  for  values  of  t outside  this  interval. 
This  is  consistent  with  the  interpretation  that  values  of  l clustered 
around  unity  arise  from  small  fluctuations  of  the  initial  perturbation 
amplitude  about  some  mean  value.  Significantly  larger  fluctuations 
cannot  be  handled  by  the  linear  model  presented  here. 


V.  CONCLUSION 

In  the  first  part  of  this  paper,  we  have  calculated  the  effect  of 
stationary  gaussian  noise  with  zero  mean  on  unstable  linear  systems. 

The  random  noise  was  considered  an  initial  condition  rather  than  a 
forcing  function.  The  development  paralleled  that  of  the  response  of  a 
linear  filter  to  random  signals.  These  results  were  specialized  to  the 
important  case  where  the  initial  disturbance  is  white  noise,  and  where 
the  instability  preferentially  amplifies  a narrow  band  of  the  Fourier 
spectrum.  This  approach  is  used  to  describe  the  statistics  of  the 
breakup  length  for  turbulent  liquid  jets  in  which  turbulence  is  viewed 
as  a random  perturbation  to  capillary  instability.  This  model  agrees 
reasonably  well  with  the  data  of  Phinney1 . 
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Figure  1.  The  statistical  distribution  of  the  normalized 
breakup  length  £ is  given  (in  arbitrary  units) 
as  a function  of  l. 
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LIST  OF  SYMBOLS 


a 

A(t) 

d 

o 

f 

F(k,t) 

6(k,t) 

i 

V Ji 

k 

L 

m,n 

N 

N 

w 

R(y,t) 

t 

u 

o 

x,y 

A 

A 

m,n 

A 

C(x,t) 

P 

o2(t) 


equilibrium  radius  of  jet 

amplitude  function,  equation  (19) 

jet  diameter;  d = 2a 
J o 

probability  density  function 

time-dependent  Fourier  transform  of  impulse  response,  Eq.  (3) 
time-dependent  power  spectrum,  Eq.  (7) 
imaginary  unit;  i = / -1 

modified  Bessel  functions  of  the  first  kind 
wave  number 

jet  breakup  length,  Eq.  (26) 
indices;  m,n  = 1,  2,  3,  ...  , N 
upper  bound 

pd  u 2 

Weber  number;  N = — = 

w T 

autocorrelation  function,  Eq.  (4) 

time  coordinate 

jet  velocity 

(same)  space  coordinate 

wavelength 

elements  of  the  covariance  matrix,  Eq.  (11) 
covariance  matri:  , Eq.  (12) 
envelope  of  ip(x,t),  Eq.  (23) 
fluid  density 

variance  of  stochastic  process;  a2(t)  = R(o,t) 
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